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A semiclassical method of complex trajectories for the calculation of the tunneling exponent in 
systems with many degrees of freedom is further developed. It is supplemented with an easily 
implementable technique, which enables one to single out the physically relevant trajectory from 
the whole set of complex classical trajectories. The method is applied to semiclassical transitions 
of a bound system through a potential barrier. We find that the properties of physically relevant 
complex trajectories are qualitatively different in the cases of potential tunneling at low energy and 
dynamical tunneling at energies exceeding the barrier height. Namely, in the case of high energies, 
the physically relevant complex trajectories describe tunneling via creation of a state close to the 
top of the barrier. The method is checked against exact solutions of the Schrodinger equation in a 
quantum mechanical system of two degrees of freedom. 
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I. INTRODUCTION 

Semiclassical methods provide a useful tool for the 
study of nonperturbative processes. Tunneling phe- 
nomena represent one of the most notable cases where 
semiclassical techniques are used to obtain otherwise 
unattainable information on the dynamics of the tran- 
sition. A standard example of semiclassical technique is 
WKB approximation to tunneling in quantum mechan- 
ics of one degree of freedom. In this case solutions S{q) 
to the Hamilton-Jacobi equation are pure imaginary in 
the classically forbidden region. Therefore, the function 
S{q) can be obtained as the action functional on a real 
trajectory ^(t), which is the solution to the equations 
of motion in Euclidean time domain, t = —it, with real 
Euclidean action Se = —iS. 

This simple picture of tunneling is no longer valid for 
systems with many degrees of freedom, where solutions 
S'(q) to the Hamilton-Jacobi equation are known to be 
generically complex in the classically forbidden region 
(see Refs. 0,13 for recent discussion). This leads to the 
concept of "mixed" tunneling, as opposed to "pure" tun- 
neling where <S'(q) is pure imaginary. "Mixed" tunnel- 
ing cannot be described by any real tunneling trajectory. 
However, it could be related to a complex trajectory. If 
so, the function S'(q) (and therefore the exponential part 
of the wave function) is calculated as the action func- 
tional on this complex trajectory. 

A particularly difficult situation arises when one con- 
siders transitions of a non-separable system with a strong 
interaction between its degrees of freedom, such that the 
quantum numbers of the system change considerably dur- 
ing the transition. Methods based on adiabatic expansion 
are not applicable in this situation, while the method of 
complex trajectories proves to be extremely useful. 

The method of complex trajectories in the form suit- 
able for the calculation of ^'-matrix elements was for- 
mulated and checked by direct numerical calculations 
in Refs. IS B, B (see Ref. j6j| for review). Further 
studies UliM C3 in 113 showed that this method 
can be generalized to the calculation of the tunneling 
wave functions and tunneling probabilities, energy split- 
tings in double well potentials and decay rates from 
metastable states. Similar methods were successful in 
the study of tunneling in high-energy collisions in field 
theory "13,1111113, where one considers systems with 
definite particle number [M — 2) in the initial state; in 
the study of chemical reactions and atom ionization pro- 
cesses, where the initial bound systems are in definite 
quantum states fe, 17, 18j; etc. The main advantage 
of the method of complex trajectories is that it can be 
easily generalized and numerically implemented in the 
cases of large and even infinite (field theory) number of 
degrees of freedom, in contrast to other methods such as 
Huygens-type construction of Refs. [u | a| and initial value 
representation (IVR) of Refs. 0, iMlfl, El E IH . 

In this paper we develop the method of complex trajec- 
tories further. Namely, we concentrate on the following 



problem. It is known [3 that the physically relevant com- 
plex trajectory satisfies the classical equations of motion 
with certain boundary conditions. However, this bound- 
ary value problem generically has also an infinite, though 
discrete, set of unphysical solutions. In one dimensional 
quantum mechanics all solutions can easily be classified. 
In systems with many degrees of freedom such a classi- 
fication is extremely difficult, if at all possible. In the 
case of small number of degrees of freedom (realistically, 
iV = 2), one can scan over all solutions and find the so- 
lution giving the largest tunneling probability 0, 0, ^3 ) 
but in systems with large or infinite number of degrees of 
freedom the problem of choosing the physically relevant 
solution becomes a formidable task. 

The problem of choosing the appropriate solution be- 
comes even more pronounced when the qualitative prop- 
erties of the relevant complex trajectory are different in 
different energy regions. This may happen when the 
physically relevant classical solution "meets" an unphys- 
ical one at some value of energy i? = E'l, or in other 
words, when solutions to the boundary value problem, 
viewed as functions of energy, bifurcate at E — Ei. 

In this paper we give an example of this sort, which 
appears to be fairly generic (see also [lU Il2. llSl lla. l2^ ) . 
We then develop a method which chooses the phyically 
relevant solution automatically, implement it numerically 
and check this method against the numerical solution to 
the full Schrodingcr equation. 

We study inelastic transitions of a bound system 
through a potential barrier. For concreteness we con- 
sider a model with one internal degree of freedom besides 
the center-of-mass coordinate. We consider a situation 
in which the spacing between the levels of the bound 
system is small compared to the height of the barrier, 
and assume strong enough coupling between the degrees 
of freedom, to make sure that the quantum numbers of 
the bound system change considerably during the transi- 
tion process. This is precisely the situation in which the 
method of complex trajectories shows its full strength. 

Transitions of bound systems involve a particular en- 
ergy scale — the height of the barrier Vq. At energies 
below Vq classical over-barrier transitions are forbidden 
energetically; the corresponding regime is called "poten- 
tial tunneling" . For £^ > Vq it is energetically allowed for 
the system to evolve classically to the other side of the 
barrier. However, over-barrier transitions may be for- 
bidden dynamically even at E > Vq. Indeed, inelastic 
interactions of a bound system with a potential barrier 
generally lead to the excitation of the internal degrees 
of freedom with the simultaneous decrease of the center- 
of-mass energy, and this may prevent the system from 
the over-barrier transition. Tunneling regime at energies 
exceeding the barrier height is called "dynamical tunnel- 
ing" 1. 



It is clear that the properties of transitions of a bound system at 
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Examples of dynamical tunneling are well-known in 
scattering theory 0. This type of tunneling between 
bound states was discovered in Ref. the generality 
of dynamical tunneling in large molecules was stressed 
in Refs. 22,, ,27J- It is dynamical tunneling that is of 
primary interest in our study. 

A novel phenomenon we observe is that dynamical tun- 
neling at E >Vo (more precisely, at E > Ei, where Ei 
is somewhat larger than Vq) occurs in the following way: 
the system jumps on top of the barrier, and restarts its 
classical evolution from the region near the top. From 
the physical viewpoint, this is not quite what is normally 
meant by "tunneling through a barrier" . Yet the transi- 
tions remain exponentially suppressed, but the reason is 
different: to jump above the barrier, the system has to 
undergo considerable rearrangement, unless the incom- 
ing state is chosen in a special way (see footnote, above). 
This rearrangement costs exponentially small probabil- 
ity factor. We note that similar exponential factor was 
argued to appear in various field theory processes with 
multi-particle final states [2II29II30II3H. 

We find that the new physical behaviour of the system 
is related to a bifurcation of the family of the complex- 
time classical solutions, viewed as functions of energy. 
This is precisely the bifurcation which we alluded to 
above. Our method of dealing with this bifurcation is 
to regularize the boundary value problem in such a way 
that the bifurcations disappear altogether (at real ener- 
gies) , and the only solutions recovered after removing the 
regularization are physical ones. 

The paper is organized as follows. The system we dis- 
cuss in this paper is introduced in Sec. Ill Al In Sec. IIIBI 
we formulate the boundary value problem for the calcu- 
lation of the tunneling exponent. Then we examine the 
classical over-barrier solutions and find all initial states 
that lead to classically allowed transitions in Sec. Ill CI 
In Sec. Ill Dl we present a straightforward application of 
the semiclassical technique, outlined in Sec. lIIIjl and find 
that it ceases to produce relevant complex trajectories n 
a certain region of initial data, namely, at E > Ei. In 
Sec. mil we introduce our regularization technique and 
show that it indeed enables one to find all the rele- 
vant complex trajectories, including ones with E > Ei 
fSec. lIIlXl . We check our method against the numerical 
solution of the full Schrodinger equation in Sec. IIII Bl In 
Sec. nrrni and Appendix O we show how our regulariza- 
tion technique is used to join smoothly the "classically 
allowed" and "classically forbidden" families of solutions 



E > Vq depend on the choice of the initial state. Namely, there 
always exists a certain class of states, transitions from which 
are not exponentially suppressed. To construct an example, one 
places the bound system on top of the barrier and evolves it 
classically backwards in time to the region where the interaction 
with the barrier is negligibly small. On the other hand, even at 
E > Vo there are states, transitions from which are exponentially 
suppressed (dynamical tunneling). 




FIG. 1: Oscillator hitting a potential barrier, with which only 
the "dark" particle interacts. 

in the cases of two- and one- dimensional quantum me- 
chanics, respectively. 



II. SEMICLASSICAL TRANSITIONS 
THROUGH A POTENTIAL BARRIER 

A. The model 

The situation we discuss in this paper is a transi- 
tion through a potential barrier of a bound system of 
Refs. ^3E3li namely a system made of two particles of 
identical mass m, moving in one dimension and bound by 
a harmonic oscillator potential of frequency cu (Fig. 
One of the particles interacts with a repulsive potential 
barrier. The potential barrier is assumed to be high and 
wide, while spacing between the oscillator levels is much 
smaller than the barrier height Vq. The Hamiltonian of 
the model is 



H 



pf pi mui , 

— + — 4 {xx 

2m 2m 4 



X2 



)2 + y^e--?/2.\ (1) 



where the conditions on the oscillator frequency and po- 
tential barrier are 



cr > h/^JmVo 



(2) 



Since the variables do not separate, this is certainly a 
non-trivial system. 

We choose units with h — 1, m — \. It is also conve- 
nient to treat the frequency w as a dimensionless param- 
eter, so that all physical quantities are dimensionless. In 
our subsequent numerical study we use the value uj — 0.5, 
still keeping, however, notation "w" in formulas. The sys- 
tem is semiclassical, i.e. conditions Q are satisfied, if one 
chooses a — \/ \/2A, Vb = 1/A, where A is a small parame- 
ter. At the classical level, this parameter is irrelevant: af- 
ter rescaling the variables^ xi xi/^/\^ xi xij-sfX^ 



^ To keep notations simple, we use the same symbols xi, X2 for 
the rescaled variables. 
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the small parameter enters only through the overall mul- 
tiplicative factor 1/A in the Hamiltonian. Therefore, the 
semiclassical technique can be developed as an asymp- 
totic expansion in A. 

The properties of the system are made clearer by re- 
placing the variables xi, X2 with the center-of-mass co- 
ordinate X = (xi + 2:2) /\/2 and the relative oscillator 
coordinate y = (xi — xn)l\f2. In terms of the latter 
variables, the Hamiltonian takes the form 



pI 



1 



-X(X+yf/2 



2 2 2 
The interaction potential 



(3) 



int 



- A 



vanishes in the asymptotic regions X — > ±00 and de- 
scribes a potential barrier between these regions. At 
X ±00 the Hamiltonian © corresponds to an os- 
cillator of frequency uj moving along the center-of-mass 
coordinate X. The oscillator asymptotic state is char- 
acterized by its excitation number N and total energy 
E = p'j^/2 + iLj{N +1/2). We are interested in the trans- 
missions through the potential barrier of the oscillator 
with given initial values of E and N. 



B. T/6 boundary value problem 

The probability of tunneling from a state with fixed 
initial energy E and oscillator excitation number N from 
the asymptotic region X — > — cx) to any state in the other 
asymptotic region X — s- +00 takes the following form: 



T{E,N) 



lim 



X 

/ 



(/|e-''-ff(*/-*0|^^7V) 



(4) 



where it is implicit that the initial and final states have 
support only well outside the range of the potential, with 
X < and X > 0, respectively. Semiclassical meth- 
ods are applicable when the initial energy and excitation 
number are parametrically large, 

E^E/X, N = N/\, 

where E and N are held constant as A ^ 0. The transi- 
tion probability has the exponential form 



T = De 



(5) 



where D is a pre-exponential factor, which is not con- 
sidered in this paper. Our purpose is to calculate the 
leading semiclassical exponent F{E, N). The expo- 
nent for tunneling from the oscillator ground state is 
obtained [H IH tllllll by taking the limit iV ^ in 
FiE,N). 

In what follows we rescale the variables, X ^ X/ \/A, 
y —> y/\/X, and omit tilde over the rescaled quantities 
E, N. 




FIG. 2; Contour in the complex time plane. 

The exponent F{E, N) is related to a complex tra- 
jectory, which satisfies a certain complexified classical 
boundary value problem. We present the derivation of 
this problem in Appendix 1X1 The outcome is as follows. 
There are two Lagrange multipliers, T and 9, which are 
related to the parameters E and N characterizing the 
incoming state. The boundary value problem is conve- 
niently formulated on the contour ABCD in the complex 
time plane (see Fig. with the imaginary part of the 
initial time equal to T/2. The coordinates X{t), y{t) 
must satisfy the complexified equations of motion in the 
internal points of the contour, and be real in the asymp- 
totic future (region D): 

, (6a) 



SX{t) 5y{t) 



Imy(i) - 
ImX(i) 



0, 

^0, 



as i — > -l-oo 



(6b) 



In the asymptotic past (region A of the contour, where 
t = t' + iT/2, t' is real negative) one can neglect the 
interaction potential C/jnt, and the oscillator decouples: 

-^(ue"''"*' +we''"*'). 



The boundary conditions in the asymptotic past, t' — > 
—00, are that the center-of-mass coordinate X must be 
real, while the complex amplitudes of the decoupled os- 
cillator must be linearly related, 



IraX ~* 0, 



as t' — c 



(6c) 



The boundary conditions l|6b|l and make, in fact, 
eight real conditions (since, e.g., lmX{t') implies 
that both ImX and ImX tend to zero), and completely 
determine a solution, up to time translation invariance 
(see discussion in Appendix 

It is shown in Appendix ^ that a solution to this 
boundary value problem is an extremum of the functional 

F[X, y- X\y*-T, 0] = -iS[X, y\ + iS[X\y*] 

- ET - N9 + Boundary terms . (7) 

The value of this functional at the extremum gives the 
exponent for the transition probability (up to large over- 
all factor 1/A, see eq. (O), 



F{E, N) = 2ImS'o(r, 



ET - NO 



(8) 
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where S'o is the action of the solution, integrated by parts, 

-^cjY~U,nt{X,y)y (9) 

Here the integration runs along the contour ABCD. The 
values of the Lagrange multipliers T and 9 are related to 
energy and excitation number as follows, 

E{T,e) = A2Im5o(r,0) , (10) 

N{T,e) = ^2lmSo{T,0) . (11) 

Making use of Eq. ||SJ|, it is straightforward to check also 
the inverse Legendrc transformation formulas, 

T{E, N) = ~-^F{E, N) , (12) 

e{E, N) = -^F{E, N) . (13) 

One can also check that the right hand side of Eq. H10() 
coincides with the energy of the classical solution, while 
the right hand side of Eq. Ijllll is equal to the classical 
counterpart of the occupation number, 

^2 

E=—+ujN; N^uv. (14) 

So, one may either search for the values of T and 9 that 
correspond to given E and N, or, following a compu- 
tationally simpler procedure, solve the boundary value 
problem © for given T and 9 and then find the corre- 
sponding values of E and N from Eq. H14|) . Note that 
the initial conditions l(Bc|) complemented by Eqs. H14|l are 
equivalent to the initial conditions of Refs. the 
latter being expressed in terms of action-angle variables. 
The boundary conditions in the asymptotic future (|6b|) 
are different from those of Refs. [3,13,13], since we consider 
inclusive, rather than fixed, final state. 

Let us discuss some subtle points of the boundary value 
problem First, one notices that the condition of 

asymptotic reality l|6b|l does not always coincide with the 
condition of reality at finite time. Of course, if the solu- 
tion approaches the asymptotic region X — > +oo on the 
part CD of the contour, the asymptotic reality condition 
(|6b|) implies that the solution is real at any finite positive 
t. Indeed, the oscillator decouples a,s X ^ +oo, so the 
condition (j6b|l means that its phase and amplitude, as 
well as X{t), are real as t — > -t-oo. Due to equations of 
motion, X{t) and y{t) are real on the entire CD-part of 
the contour. This situation corresponds to the transition 
directly to the asymptotic region X +oo. However, 
the situation can be drastically different if the solution 
on the final part of the time contour remains in the inter- 
action region. For example, let us imagine that the solu- 
tion approaches the saddle point of the potential X = 0, 



y = a,s t ^ +00. Since one of the perturbations about 
this point is unstable, there may exist solutions which 
approach this point exponentially along the unstable di- 
rection, i.e. X{t), y{t) oc e~'^°"''* * with possibly complex 
pre- factors. In this case the solution may be complex at 
any finite time, and become real only asymptotically, as 
t — > +00. Such solution corresponds to tunneling to the 
saddle point of the barrier, after which the system rolls 
down classically towards X +oo (with probability of 
order 1, inessential for the tunneling exponent F). We 
will see in Sec. IIII Al that the situation of this sort in- 
deed takes place for some values of energy and excitation 
number. 

Second, since at large negative time (in the asymptotic 
region X — > —oo) the interaction potential disappears, 
it is straightforward to continue the asymptotics of the 
solution to the real time axis. For solutions satisfying 
l|Hc|l this gives at large negative time 

y{t) - ( ue--^/2g-»^* _^ ^*^9+uT/2^ru.t\ 

lmX{t)^-^Px . 

We see that the dynamical coordinates on the negative 
side of the real time axis are generally complex. For 
solutions approaching the asymptotic region X +oo 
as t +O0 (so that X and y are exactly real at finite 
t > 0), this means, that there should exist a branch point 
in the complex time plane: the contour A'ABC in Fig. [21 
winds around this point and cannot be deformed to the 
real time axis. This argument does not work for solutions 
ending in the interaction region as t -\-oo, so branch 
points between the AB-part of the contour and the real 
time axis may be absent. We will see in Sec. IIII Al that 
this is indeed the case in our model in a certain range of 
E and N. 



C. Over-barrier transitions: the region of 
classically allowed transitions and its boundary 

Eo{N) 

Before studying the exponentially suppressed transi- 
tions, let us consider the classically allowed ones. To 
this end, let us study the classical evolution (real time, 
real- valued coordinates), in which the system is initially 
located at large negative X and moves with positive 
center-of-mass velocity towards the asymptotic region 
X +O0. The classical dynamics of the system is speci- 
fied by four initial parameters. One of them (e.g., the ini- 
tial center-of-mass coordinate) fixes the invariance under 
time translation, while the other three are the total en- 
ergy E, the initial excitation number of the y-oscillator, 
defined in classical theory as = Eqsc/^i and the initial 
oscillator phase ipi. 

Any initial quantum state of our system can be fully 
determined by energy E and initial oscillator excitation 
number N] we can represent each state by a point in the 
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E-N plane. There is, however, one additional classically 
relevant initial parameter, the oscillator phase (pi. An 
initial state {E, N) leads to unsuppressed transmission if 
the corresponding classical over-barrier transitions^ are 
possible for some value(s) of ipi. These states form some 
region in the E-N plane, which is to be found in this 
section. 

For given iV, at large enough E the system can cer- 
tainly evolve to the other side of the barrier. On the other 
hand, if E is smaller than the barrier height, the system 
definitely undergoes reflection. Thus, there exists some 
boundary energy Eq{N), such that classical transitions 
are possible for E > Eo{N), while for E < EoiN) they do 
not occur for any initial phase ipi. The line Eq(N) rep- 
resents the boundary of the region of classically allowed 
transitions. We have calculated Eo{N) numerically: the 
result** is shown in Fig. 13 

An important point of the boundary Eo{N) corre- 
sponds to the static unstable classical solution X{t) — 
y{t) = 0. In field theory context, such a solution is 
called "sphaleron" [s^]), and for the sake of terminol- 
ogy we will use this name in what follows. This so- 
lution is the saddle point of the potential U{X,y) = 
uj'^y^/2 + Uint{X,y) and has exactly one unstable direc- 
tion, the negative mode (see Fig. EJ. The sphaleron en- 
ergy Eg = [/(O, 0) = 1 determines the minimum value of 
the function Eo{N). Indeed, while classical over-barrier 
transitions with E < Es are impossible, the over-barrier 
solution with slightly higher energy can be obtained as 
follows: at the point X — y — 0, one adds momentum 
along the negative mode, thus "pushing" the system to- 
wards X — ^ -|-oo. Continuing this solution backwards 
in time one finds that the system tends to X — > — oo 
for large negative time and has a certain oscillator ex- 
citation number. Solutions with energy closer to the 
sphaleron energy correspond to smaller "push" and thus 
spend longer time near the sphaleron. In the limiting 
case when the energy is equal to Eg , the solution spends 
an infinite time in the vicinity of the sphaleron. This 
limiting case has a definite initial excitation number iVs, 
so that Eo{Ns) = Es (see Fig. EJ. The value of iVs is 
unique because there is exactly one negative direction of 
the potential in the vicinity of the sphaleron. 

In complete analogy to the features of the over-barrier 
classical solutions near the sphaleron point {Eg, Ng), one 
expects that as the values of E, N approach any other 
boundary point {Eo{N), N), the corresponding over- 



^ Note that the corresponding classical solutions obey the bound- 
ary conditions 16 bi . I6cl with T = 9 = 0, i.e., they are solutions 
to the boundary value problem jSJ. 

* Note that the boundary Eo{N) of the region of classically allowed 
transitions can be extended to Af > A^g. As E = Eg is the 
absolute minimum of the energy of classically alowed transitions, 
the function Eo{N) grows with N at N > Nq. In fact, it tends 
to asymptotics Eq^ = ujN as Af — > -|-oo. In what follows we are 
not interested in transitions with A'^ > A^g, so this part of the 
boundary Eo{N) is not presented in Fig. 1^ 




0.5 Eg = 1 1.5 

E 



FIG. 3: The boundary Eo(N) of the region of classically al- 
lowed transitions, the bifurcation line Ei{N) and the line of 
the periodic instantons Epi{N). 

barrier solutions will spend more and more time in the 
interaction region, where C/int ^ 0. This follows from 
a continuity argument. Namely, let us first fix the ini- 
tial and final times, and tf. If in this time interval 
a solution with energy Ei evolves to the other side of 
the barrier and a solution with energy E2 and the same 
oscillator excitation number is reflected back, there ex- 
ists an intermediate energy at which the solution ends 
up at t = in the interaction region. Taking the limit 
— > -foo and {Ei — i?2) — 0, we obtain a point at the 
boundary Eo{N) and a solution tending asymptotically 
to some unstable time dependent solution that spends 
infinite time in the interaction region. We call the latter 
solution excited sphaleron; it describes some (in general 
nonlinear) oscillations above the sphaleron along the sta- 
ble direction in coordinate space. Therefore, every point 
of the boundary {Eo{N), N) corresponds to some ex- 
cited sphaleron. Solutions tending asymptotically to the 
excited sphalerons, form a surface in a phase space (sepa- 
ratrix), which separates regions of qualitatively different 
classical motion of the system. 

We display in Fig. ^ a solution, found numerically in 
our model, that tends to an excited sphaleron. We see 
that the trajectory of excited sphaleron is, roughly speak- 
ing, orthogonal to the unstable direction at the saddle 
point {X = 0,y = 0). 

D. Suppressed transitions: bifurcation line Ei{N) 

Let us now turn to classically forbidden transitions, 
and consider the boundary value problem ©. It is rela- 
tively straightforward to obtain numerically solutions for 
6 = 0. The boundary conditions Ij6bp . (|Sc|l in this case 
take the form of reality conditions in the asymptotic fu- 
ture and past. It can be shown ;3 that the physically 
relevant solutions with 6 = are real on the entire con- 
tour ABCD of Fig. 121 and describe nonlinear oscillations 
in the upside-down potential on the Euclidean part BC 
of the contour. The period of the oscillations is equal to 
T, so that the points B and C are two different turning 
points, where X = y = 0. These real Euclidean solu- 
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y-X 



y+X 



FIG. 4: The potential (dotted lines) in the vicinity of the 
sphaleron {X = 0, y ~ 0) (marked by the point), excited 
sphaleron (thick line) corresponding to the point {E, N) — 
(1.985,3.72). at the boundary of the region of classically al- 
lowed transitions, and the trajectory of the solution which is 
close to this excited sphaleron (thin line). In this figure the 
asymptotic regions X —> ±cxd are along the diagonal. 



tions are called periodic instantons. A practical tech- 
nique for obtaining these solutions numerically on the 
Euclidean part BC consists in minimizing the Euclidean 
action (for exa mple, w ith the method of conjugate gra- 
dients, see Ref. [ill [12 for details). The solutions on the 
entire contour are then obtained by solving numerically 
the Cauchy problem, forward in time along the line CD 
and backward in time along the line BA. Having the so- 
lution in asymptotic past (region A), one then calculates 
its energy and excitation number H14I) . The solutions to 
this Cauchy problem are obviously real, so the bound- 
ary conditions Ij6b|l . (|Sc|l are indeed satisfied with 6 = 0. 
It is worth noting that solutions with 9 = are similar 
to the ones in quantum mechanics of one degree of free- 
dom. The line of periodic instantons in E—N plane in 
our model is shown in Fig. |31 

Once the solutions with 6 = are found, it is natural 
to try to cover the entire region of classically forbidden 
transitions oi E~N plane with a deformation procedure, 
by moving in small steps in 6 and T . The solution to the 
boundary value problem with (T -|- AT, 6 + lS.6) may be 
obtained numerically, by applying an iteration technique, 
with the known solution at (T, 6) serving as the initial 
approximation^. Provided the solutions end up in correct 
asymptotic region at each step, i.e. X — > -\-oo on part D 
of the contour, the solutions obtained by this procedure 



10 15 20 25 



^ In practice, the Newton-R aph son method is particularly conve- 
nient (see Refs. [TllTllTllllV 



FIG. 5: The dependence of the tunneling coordinate X on 
time for two solutions with nearly the same energy and ini- 
tial excitation number. The physical solution tunnels to 
the asymptotic region X ^ -|-oo, while the unphysical one 
gets reflected back to X —* ~oo. The physical solution 
has E = 1.028, A*' = 0.44, while the unphysical one has 
E = 1.034, A'' = 0.44. These two solutions are close to the 
point at the bifurcation line Ei{N = 0.44) = 1.031. 



of small deformations are physically relevant. However, 
the method of small deformations fails to produce rele- 
vant solution if there are bifurcation points in the E-N 
plane, where the physical branch of solutions merges to 
an unphysical branch. As there are unphysical solutions 
close to physical ones in the vicinity of bifurcation points, 
one cannot use the procedure of small deformations near 
these points. 

We have found numerically that in our model the 
method of small deformations produces correct solutions 
to the T/6 boundary value problem in a large region of 
the E — N plane, where E < Ei{N). However, at sufh- 
ciently high energy E > Ei{N), where Ei{N) > Es, the 
deformation procedure generates solutions that bounce 
back from the barrier (see Fig.O, i.e. have wrong "topol- 
ogy" . This happens deep inside the region of classically 
forbidden transitions, where the suppression is large, and 
one naively expects the semiclassical technique to work 
well. Clearly, solutions with wrong topology do not de- 
scribe the tunneling transitions of interest. Therefore, if 
the semiclassical method is applicable at all in the re- 
gion Ei{N) < E < Eo{N), there exists another, physical 
branch of solutions. In that case the line Ei{N) is the 
bifurcation line, where the physical solutions "meet" the 
ones with wrong "topology" . Walking in small steps in 
6 and T is useless in the vicinity of this bifurcation line, 
and one needs to introduce some trick to find the relevant 
solutions beyond that line. The bifurcation line Ei{N) 
for our quantum mechanical problem of two degrees of 
freedom is shown in Fig. O 

The loss of topology beyond a certain bifurcation line 
m E — N plane is by no means a property of our model 
only. This phenomenon has been observed in field theory 
models, in the context of both induced false vacuum de- 
cay jl4j | and baryon-number violating transitions in gauge 
theory p5j (in field-theoretic models, the parameter N is 
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the number of incoming particles). In all cases, the loss 
of topology prevented one to compute the semiclassical 
exponent for the transition probability in an interesting 
region of relatively high energies. 

Coming back to quantum mechanics of two degrees of 
freedom, we point out that the properties of tunneling 
solutions with different energies approaching the bifurca- 
tion line Ei{N) from the left of the E-N plane, are in 
some sense similar to the properties of tunneling solutions 
in one-dimensional quantum mechanics whose energy is 
close to the barrier height, see Appendix |0 Again by 
continuity, these solutions of our two-dimensional model 
spend a long time in the interaction region; this time 
tends to infinity on the line Ei(N). Hence, at any point 
of this line, there is a solution which starts in the asymp- 
totic region left of the barrier, and ends up on an excited 
sphaleron. Such a behavior is indeed possible because of 
the existence of an unstable direction near the (excited) 
sphaleron, even for complex initial data. We suggest in 
the next Section a trick to deal with this situation — this 
is our regularization technique. 



III. REGULARIZATION TECHNIQUE 

In this Section we develop our regularization technique, 
and find the physically relevant solutions between the 
lines Ei{N) and Eo{N). We will see that all solutions 
from the new branch (and not only on the lines Eo{N) 
and Ei{N)) correspond to tunneling onto the excited 
sphaleron ("tunneling on top of the barrier"). These so- 
lutions would be very difficult, if at all possible, to ob- 
tain directly, by solving numerically the non-regularized 
classical boundary value problem ©: they are complex 
at finite times, and become real only asymptotically as 
t -|-oo, whereas numerical methods require working 
with finite time intervals. 

As an additional advantage, our regularization tech- 
nique enables one to obtain a family of the over-barrier 
solutions, that covers all the region of the initial data, 
corresponding to classicaly allowed transitions, including 
its boundary. This is of interest in models with large 
number of degrees of freedom and in field theory, where 
finding the boundary Eq (N) by direct methods is difficult 
(see e.g., Ref. for discussion of this point). 



A. Regularized problem: classically forbidden 
transitions 

The main idea of our method is to regularize the equa- 
tions of motion by adding a term proportional to a small 
parameter e so that configurations staying for an infi- 
nite time near the sphaleron no longer exist among the 
solutions of the T/9 boundary value problem. After per- 
forming the regularization we explore all region of classi- 
cally forbidden transitions without crossing the bifurca- 
tion line. Taking then the limit e ^ we reconstruct the 



correct values of F, E and N. 

When formulating the regularization technique it 
is more convenient to work with the functional 
F[X,y;X*,y*;T,9], Eq. Q, itself rather than with the 
equations of motion. We prevent F from beeing extrem- 
ized by configurations approaching the excited sphalerons 
asymptotically. To achieve this, we add to the original 
functional Q a new term of the form 2eTiat^ where Tint 
estimates the time the solution "spends" in the inter- 
action region. The parameter of regularization e is the 
smallest one in the problem, so any regular extremum 
of the functional F (the solution that spends finite time 
in the region Uint ^ 0) changes only slightly after the 
regularization. At the same time, the excited sphaleron 
configuration has Tint = oo which leads to infinite value 
of the regularized functional F^ = F + 2eTint- Hence, 
the excited sphalerons are not stationary points of the 
regularized functional. 

For the problem at hand ?7int ^ 1 in the interaction 
region, and Tint can be defined as follows. 



T 

Jint — 77 



dtUint{X,y)+ / dtUi,,t{X*,y*) 



(15) 



We notice that Tint is real, and that the regularization is 
equivalent to the multiplication of the interaction poten- 
tial by a complex factor 



(1 - «e)f/ir 



(16) 



This results in the corresponding change of the clas- 
sical equations of motion, while the boundary condi- 
tions l|6b|l . (|Sc|l remain unaltered. 

We still have to understand whether solutions with 
e 7^ exist at all. The reason for the existence of such 
solutions is as follows. Let us consider a well-defined (for 
e > 0) matrix element 



lim 



where \E, N) denotes, as before, the incoming state with 
given energy and oscillator excitation number. The quan- 
tity 7^ has a well defined limit as e ^ 0, equal to the 
tunneling probability 10} . As the saddle point of the reg- 
ularized functional T^ = F + 2eTint gives the semiclassical 
exponent for the quantity 7^ , one expects that such sad- 
dle point indeed exists. 

Therefore, the regularized T /6 boundary value prob- 
lem is expected to have solutions necessarily spending 
finite time in the interaction region. By continuity, these 
solutions do not experience reflection from the barrier, if 
one makes use of the procedure of small deformations 
starting from solutions with correct "topology" . The 
line Ei{N) is no longer a bifurcation line of the regu- 
larized system, so the procedure of small deformations 
enables one to cover the entire region of classically for- 
bidden transitions. The semiclassical suppression factor 
of the original problem is recovered in the limit e — > 0. 
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It is worth noting that the interaction time is Legendre 
conjugate to e, 



ld_ 

2d^ 



Fe{E,N,e) 



(17) 



This equation may be used as a check of numerical cal- 
culations. 

We implemented the regularization procedure numeri- 
cally. To solve the boundary value problem, we make use 
of the computational methods described in Ref. 0, 0| . 
To obtain the semiclassical tunneling exponent in the re- 
gion between the bifurcation line Ei (N) and the bound- 
ary of the region of classically allowed transitions Eq (N) , 
we began with a solution to the non-regularized problem 
deep in the "forbidden" region of initial data (i.e., at 
E < Ei{N)). Then we increased the value of e from zero 
to a certain small positive number, keeping T and fixed. 
Then we changed T and 6 in small steps, keeping e finite, 
and found solutions to the regularized problem in the re- 
gion Ei{N) < E < Eq{N). These solutions had correct 
"topology" , i.e. they indeed ended up in the asymptotic 
region X — > -|-oo. Finally, we lowered e and extrapolated 
F, E and N to the hmit e ^ 0. 

Let us consider more carefully the solutions in the re- 
gion Ei{N) < E < Eo{N) which we obtain in the limit 
e — > 0. They belong to a new branch, and thus may 
exhibit new physical properties. Indeed, we found that, 
as the value of e decreases to zero, the solution at any 
point {E,N) with Ei{N) < E < Eo{N) spends more 
and more time in the interaction region. The limiting 
solution corresponding to e = has infinite interaction 
time: in other words, it tends, as < — > +oo, to one of 
the excited sphalerons. The resulting physical picture is 
that at large enough energy (i.e., at > Ei{N)), the 
system prefers to tunnel exactly onto an unstable classi- 
cal solution, excited sphaleron, that oscillates about the 
top of the potential barrier. To demonstrate this, we 
have plotted in Fig. the solution x{t) = {X{t),y{t)) 
at large times, after taking numerically the limit e — > 0. 
To understand this figure, one recalls that the potential 
near the sphaleron point X = y = has one positive 
mode and one negative mode. Namely, by introducing 
new coordinates c+, c_, 

X = cos a c_|_ + sin a c_ , 
y = — sin a c+ + cos a c_ , 

ctg2a = — - , 
one writes, in the vicinity of the sphaleron, 

^ - ^ + 2 + 2 + 2 2 ' 

where 




±{-l 



-) + \ 1 + 



> 



FIG. 6; The large-time behavior of a solution with e = 
at {E = 1.05, N = 0.43). The coordinates X and y are 
decomposed in the basis of the eigenmodes near the sphaleron. 
Note that Im c+ = 0. 



Since the solutions to the T/9 boundary value problem 
are complex, the coordinates c+ and c_ are complex too. 
We show in Fig. |S1 real and imaginary parts of c+ and 
c_ at large real time t (part CD of the contour). We see 
that while Rec+ oscillates, the unstable coordinate c_ 
approaches asymptotically the sphaleron value: c_ — > 
as f — > -|-oo. The imaginary part of c_ is non-zero at any 
finite time. This is the reason for the failure of straight- 
forward numerical methods in the region E > Ei{N): 
the solutions from the physical branch do not satisfy the 
conditions of reality at any large but finite final time. We 
have pointed out in Sec. lIIBl that this can happen only if 
the solution ends up near the sphaleron, which has a neg- 
ative mode. This is precisely what happens: for e = at 
asymptotically large t our solutions are real and oscillate 
near the sphaleron, remaining in the interacton region. 



B. Regularization technique versus exact 
quantum— mechanical solution 

Quantum mechanics of two degrees of freedom is a 
convenient testing ground for checking the semiclassi- 
cal methods and, in particular, our regularization tech- 
nique. We have found the solutions to the full stationary 
Schrodinger equation and exact tunneling probability T 
by applying the numerical technique of Refs. 0, fl^ . 
Our numerical calculations were performed for several 
small values of the semiclassical parameter A, namely, 
for A = 0.01 — 0.1. Transitions through the barrier for 
these values of the semiclassical parameter are well sup- 
pressed. In particular, for A = 0.02 the tunneling prob- 
ability T is of order e^^**. To check the semiclassical 
result with better precision, we have calculated the exact 
suppression exponent FqM(A) = — AlogT (cf. ^) for 
A = 0.09, 0.05, 0.03, 0.02 and extrapolated Fqm to A = 
by polynomials of the third and fourth order. The re- 
sults of extrapolation are independent of the order (3 or 
4) of polynomials with precision of 1%. The extrapolated 
suppression exponent -F'qm(O) corresponds to infinite sup- 
pression and must exactly coincide (up to numerical er- 
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FIG. 7: The tunneling exponent F{E, N) in the region E > 
Es = 1- The lines show the semiclassical results while points 
represent exact ones, obtained by solving the Schrodinger 
equation. The lines across the plot are the boundary of the 
region of classically allowed transitions Eo{N) and the bifur- 
cation line Ei{N). 



tots) with the correct semiclassical result. 

We performed this check in the region E > Es — 1, 
which is most interesting for our purposes. The results 
of the full quantum mechanical calculation of the sup- 
pression exponent -Fqm in the limit A ^ are repre- 
sented by points in Fig. 13 The lines in that figure rep- 
resent the values of the semiclassical exponent F{E, N) 
for constant N, which we obtain in the limit e — + of the 
rcgularization procedure. In practice, instead of taking 
the limit e ^ one calculates the regularized functional 
F,{E,N) = F{E,N) + 0(e) for small enough e. In our 
calculations we used the value e = 10~^, so that the value 
of the suppression exponent was found with precision of 
order 10~^. We see that in the entire region of classically 
forbidden transitions (including the region E > Ei{N)) 
the senriclassical result for F coincides with the exact 
one. 



C. Classically allowed transitions 

Let us show that our rcgularization procedure enables 
one to obtain a subset of classical over-barrier solutions, 
existing at high enough energies. This subset is inter- 
esting, as it extends all the way to the boundary of the 
region of classically allowed transitions, E ~ Eo{N). In 
principle, finding this boundary is purely a problem of 
classical mechanics, and, indeed, in mechanics of two de- 
grees of freedom one obtains this boundary numerically 
by solving the Cauchy problem for given E and N and 
all possible oscillator phases, see Sec. Ill CI However, if 
the number of degrees of freedom is much larger, this 
classical problem becomes quite complicated, as one has 




(c) 

FIG. 8: The phase of the tunneling coordinate in complex 
time plane at the three points of the curve r = 380, i9 = 130. 
Points (a), (b), (c) correspond to e = ea = 0.01, e — et = 
0.0048 and e — ec = 0, respectively. The asymptotics X 
— CO and X +oo correspond to arg(X) = tt and 0. The 
contour in the time plane is plotted with white line. 



to span a high-dimensional space of Cauchy data. As an 
example, a stochastic Monte Carlo technique was devel- 
oped in Ref. [s^ to deal with this problem in field theory 
context. The approach below is as an alternative to the 
Cauchy methods. 

First, let us recall that all classical over-barrier solu- 
tions with given energy and excitation number satisfy the 
T/6 boundary value problem with T = 0, 9 = 0. We can- 
not reach the "allowed" region oi E — N plane without 
regularization, since we have to cross the line Eo{N) cor- 
responding to the excited sphaleron configurations in the 
final state. However, the excited sphalerons are longer 
exist among the solutions of the regularized boundary 
value problem at any finite value of e. This suggests that 
the regularization enables one to enter the region of the 
classically allowed transitions and, after taking an appro- 
priate limit, obtain classical solutions with finite values 
of -B, N. 

By definition, the classically allowed transitions have 
F = 0. Thus, one expects that in the "allowed" region 
of initial data, the regularized problem has the property 
that F^{E,N) = ef{E,N) + 0{e^). In view of the inverse 
Legendre formulas (|12|) , (|13|) the values of T and 6 must 
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be of order e: T = eT{E,N), 9 = ed{E,N), where the 
quantities r and are related to the initial energy and 
excitation number (see Eqs. p2|l . 113|l ') in the following 
way, 

where we have used Eq. (|17|l . Therefore, one expects that 
one can invade the region of classically allowed transi- 
tions by taking a fairly sophisticated limit e ^ with 
r = T/e = const, -d = 9/e = const. For the allowed tran- 
sitions the parameters t and ■)? are analogous to T and 

e. 

By solving the regularized T/9 boundary value prob- 
lem one constructs a single solution for given E and N. 
On the other hand, for e = there are more classical 
over-barrier solutions: they form a continuous family la- 
beled by the initial oscillator phase. Thus, after taking 
the limit e ^ one obtains a subset of over-barrier solu- 
tions, which should therefore obey some additional con- 
straint. It is almost obvious, that this constraint is that 
the interaction time Tint, Eq. I|15|) . is minimal. This is 
shown in Appendix IbI 

The subset of classical over-barrier solutions obtained 
in the limit e — > of the regularized T/9 procedure ex- 
tends all the way to the boundary of the region of classi- 
cally allowed transitions. Let us see what happens as one 
approaches this boundary from the "classically allowed" 
side. At the boundary Eo{N), the unregularized solu- 
tions tend to excited sphalerons, so the interaction time 
Tint is infinite. This is consistent with (|18l) . I|19() only if r 
and I? become infinite at the boundary. Thus, to obtain 
a point of the boundary one takes the further limit, 

{Eo{N),N)^ hm {E{t,^),N{t,^)) . 

T const 

r — *-t-oo 

Different values of t/t!} correspond to different points of 
the line Eo{N). In this way one finds the boundary of 
the region of classically allowed transitions without an 
initial-state simulation. 

We have checked this procedure numerically. The limit 
e ^ exists indeed — the values of E and N tend to 
the point of the E-N plane, which corresponds to the 
classically allowed transition. The phase of the tunneling 
coordinate X(t) in complex time plane is shown in Fig.|Sl 
for the three points (a), (b) and (c) of the curve t = 
T/e = 380, ■d = 9/€= 130. Point (a) hes deep inside the 
tunneling region, Ea < Ei{Na), point (c) corresponds to 
over-barrier solution with T = 0, 9 = 0, e = 0, point (b) 
is in the middle of the curve. The branch points of the 
solution®, the cuts and the contour are clearly seen on 
these graphs. 




E 



(b) 

FIG. 9: Dependence of the parameter T = —dF/dE on en- 
ergy for (a) two-dimensional model with fixed A'' = 0.1 and 
(b) one-dimensional model (see AppendixjUJ. Different lines 
correspond to diflferent branches of classical solutions to T/9 
boundary value problem. The branches labelled "reflection" 
end up on the wrong side of the barrier. Figure (b) contains 
also a line with nonzero e. 



It is worth noting that the left branch points move 
down as T and 9 approach zero. Solutions close enough 
to the boundary Eq{N) have left branch point in the 
lower complex half-plane, see Fig. |S1 Therefore, the cor- 
responding contour may be continuously deformed to the 
real time axis. These solutions still satisfy the reality con- 
ditions asymptotically (see Fig. El, but show nontrivial 
complex behavior at any finite time. 

Making use of the regularized T/9 procedure, one is 
able to approach the boundary of the region of classi- 
cally allowed transitions from both sides. The points 
at this boundary are obtained by taking the limits 
T ^ 0, T/9 = const of the tunneling solutions and 
T -t-oo, t/i3 = const of the classically allowed ones. 



The phase of the tunneUng coordinate turns by n around the 



branch point. The points where the phase of the tunneling coor- 
dinate turns by 27r correspond to the zeroes of X{t). 
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As T* = r/d = T/9 hy construction, the lines r* = const 
are continuous at the boundary Eq(N), though may have 
discontinuity of the derivatives. The variable r* can be 
used to parametrize the curve Eq(N). 



IV. CONCLUSIONS 

We conclude that classical solutions describing trans- 
missions of a bound system through a potential barrier 
with different values of energy and initial oscillator ex- 
citation number form three branches. These branches 
merge at bifurcation lines Eo{N) and Ei{N). Solu- 
tions from different branches describe physically differ- 
ent transition processes. Namely, solutions at low en- 
ergies E < Ei{N) describe conventional potential-like 
tunneling, while at i? > Eo{N) they correspond to un- 
suppressed over-barrier transitions. At intermediate en- 
ergies, Ei{N) < E < Eq{N), physically relevant solu- 
tions describe transitions on top of the barrier. This 
branch structure is shown in Fig. |^, where the period 
T = OF/ dE obtained numerically for the solutions from 
different branches is plotted as function of energy for 
N = 0.1. 

One notices that the qualitative structure of branches 
in model, with internal degrees of freedom is similar to 
the structure of branches in one-dimensional quantum 
mechanics (see Appendix 0. The latter is shown in 
Fig. O). The features of solutions in both cases are simi- 
lar, though the solutions ending up on top of the barrier 
are degenerate in energy in one-dimensional case, and 
hence are not really physically interesting. 

In this paper we introduced the regularization tech- 
nique which enables one to connect smoothly solutions 
from different branches. Its advantage is that it auto- 
matically chooses the physically relevant branch. This 
technique is particularly convenient in numerical stud- 
ies: we have seen that it enables one to cover the whole 
interesting region of parameter space. We applied this 
technique to baryon number violating processes in elec- 
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APPENDIX A: T/O BOUNDARY VALUE 
PROBLEM 



The semiclassical method for calculating the proba- 
bility of tunneling from a state with a few parameters 
fixed was developed in 0, F^J l Hpp i n context of field 
theoretical models and in [sl \A fsl IllL in quantum 
mechanics. Here we outline the method, adapted to our 
model of two degrees of freedom. 



1. Path integral representation of the transition 
probability 

We begin with the path integral representation for 
the probability of tunneling from the asymptotic region 
X — > —oo through a potential barrier. Let the incoming 
state \E, N) have fixed energy and oscillator excitation 
number, and has support only for X <ti 0, well outside 
the range of the potential barrier. The inclusive tunnel- 
ing probability for states of this type is given by 



dyf 



(X/,y/|e-^(*/-*')|i?,7V) , (Al) 



where H is the Hamiltonian operator. This probability 
can be reexpressed in terms of the transition amplitudes 

Af.^iXf, y/|e-'^(*^-*')|X„ y,) (A2) 
and initial-state matrix elements 

Bw = {X,, y,\E, N){E, N\Xi y[) (A3) 
in the following way, 

{+00 
j dXf J dX, dxl 
n — rv-> 



+ 00 



dyi dy[ dyf AfiA*, jB, 



(A4) 



The transition amplitude and its complex conjugate have 
the familiar path integral representation: 



= j [dx\ 
A*,f = [ [dx' 



JS[x] 



(A5) 



iSlx'' 
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where x = {X, y) , and S is the action of the model. To 
obtain a similar representation for the initial-state matrix 
elements, let us rewrite Bui as follows, 



(A6) 



where P/v and Pe denote the projectors onto states with 
oscillator excitation number N and total energy E re- 
spectively. It is convenient to use the coherent state for- 
malism for the y-oscillator and choose the momentum 
basis for the X-coordinate. In this representation, the 
kernel of the projector operator PePn takes the form 

{q, b\PEPN\p, a) = j^-:^ J d^dri 

exp (^-iE^ - iNr] + ^p^^ + e^^^+'^ba^ 5{q ~ p) , 



where 

F[x, x'- T, e] = -iS[X,v]+iS[X', y'] 

- ET - NO + B,{xi,x'f,T,e). (A9) 
Here the non-trivial initial term Bi is 



B, = 



2T 



'i(y,2 + yf)(l+e2™) 



(AlO) 



In (|A8I) x and x' are independent integration variables, 
while x'j: = Xf, see Eqs. (jA4ll . 



where \p, a) is the eigenstate of the center-of-mass mo- 
mentum px and y-oscillator annihilation operator a with 
eigenvalues p and a respectively. It is straightforward to 
express this matrix element in the coordinate represen- 
tation using the formulas 



{X\p) = 



1 



Evaluating the Gaussian integrals over a, 5, p, q, we ob- 
tain 

Bw = Jd^df] exp (^-lE^ - iN-q - li^i—^ 



+ 



(1 +e-2«'^«-2«»,) 



(A7) 



where we omit the pre-exponential factor depending on 
ry, ^. For the subsequent formulation of the boundary 
value problem it is convenient to introduce the notations 



T = -i^ 



= — 1?7 



Then, combining together the integral representations 
(|A7I) and (|A5p and rescaling coordinates, energy and ex- 
citation number a; — > x ■/%/A, E E/X, N N/\, we 
finally obtain 



J dT de j [dxdx] 



exp^--F[f, x'- T,e]\\ , (A8) 



2. The boundary value problem 



For small A, the path integral ljA8|l is dominated by 
a stationary point of the functional F. Thus, to cal- 
culate the tunneling probability exponent, we extremize 
this functional with respect to all variables of integration: 
X{t), y{t), X'{t), y'{t), T, 9. Note that because of the 
limit tf — ti — > +00, the variation with respect to the 
initial and final values of coordinates leads to boundary 
conditions imposed at asymptotic t ±oo, rather than 
at finite times ti,tf. Note also that the stationary points 
may be complex. 

The variation of the functional HA9|I with respect to 
the coordinates at intermediate times gives second order 
equations of motion, in general complexified. 



SS 



ss 



SS' 



SS' 



5X{t) 5y{t) 5X'{t) 6y'{t) 







(Alia) 



The boundary conditions at the final time — > +oo are 
obtained by extremization of F with respect to Xf = X'^ 

Vf 



y'j. These are 



Xf 



yf = Vf 



(Allb) 



It is convenient to write the conditions at the initial 
time (obtained by varying Xi, yi, X[, y'^) in terms of 
the asymptotic quantities. At the initial moment of time 
ti — > — oo, the system moves in the region X — > — oo, well 
outside of the range of the potential barrier. Equations 
(|Allap in this region describe free motion of decoupled 
oscillator, and the general solution takes the following 
form, 



X{t) = X,+p,{t~t,), 



yit) 



1 



^—iLj{t — ti) 



ae 



while the solution for X'{t), y'{t) has similar form. For 
the moment, a and a are independent variables. The 
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initial boundary conditions in terms of the asymptotic 
variables Xi, pi, a, d take the form: 



Pt = Pi 



iT ' 



(Allc) 



a + d^ a'e"'^^" 



The variation with respect to the Lagrange multipliers T 
and 6 gives the relation between the values of E, N and 
initial asymptotic variables (here we use the boundary 
conditions IjAllcjl ). 



E 

2 

N ^ ad . 



LUN 



(Alld) 



Equations (|Alla|) - IjAlldp constitute the complete set 
of saddle-point equations for the functional F. 

The variables X' and y' originate from the conju- 
gate amplitude A*,j: (see Eq. (|A5ll '). which suggests that 
they are the complex conjugate to X, y. Indeed, the 
Ansatz X'{t) — X*{t), y'{t) — y*{t) is compatible with 
the boundary value problem IjAlip . Then the Lagrange 
multipliers T, 9 are real, and the problem ljAll|l may 
be conveniently formulated at the contour ABCD in the 
complex time plane (see Fig.|2}. 

Now we have only two independent complex variables 
X{t) and y{t), that have to satisfy the classical equations 
of motion in the interior of the contour. 



6S 



SS 







(Al2a) 



SX{t) 5y{t) 

The final boundary conditions (sec Eq. ljAllb|) ) become 
the conditions of the reality of the variables X{t) and 
y{t) at the asymptotic part D of the contour: 



t 



(A12b) 



ImX/ = 0, Imy/ = , 

ImX/ = l-mijf = , 

The seemingly complicated initial conditions IjAllcfl sim- 
plify when written in terms of the time coordinate t' — 
t+iT/2 running along the part AB of the contour. Let us 
again write the asymptotics of a solution, but now along 
the initial part AB of the contour: 

X = Xo + po{t' ~ U) , 



1 



y 



ue 



-iuj{t'-ti) _|_ ^gji^(t'-ti) 



In terms of Xq, y^, u and u, the boundary conditions 
(|Allc|) are 



ImXo = 0, Impo = , 



(A12c) 



* 9 

u e 



Finally, we write Eqs. IjAlldp in terms of the asymptotic 
variables along the initial part of the contour: 



E = 

TV : 



2 

: UJUV 



-ljN 



(A13) 



These equations determine the Lagrange multipliers T, 9 
in terms of E, N . Alternatively, we can solve the prob- 
lem (IA12|I for given values of T, 9 and find the values of 
E,N from Eqs. (|IT3|) . what is more convenient compu- 
tationally. 

Given a solution to the problem ljA12p . the exponent 
F is the value of the functional (jA9|l at this saddle point. 
In this way we obtain the expression ||HJ| for the tunneling 
exponent. The exponent F is expressed now in terms of 
So, eq. Q — integrated by parts action of the system. 
Non-trivial boundary term Bi, eq. (|A10p . is cancelled by 
the boundary term coming from the integration by parts. 
Note that we did not make use of the constraints (|A13p 
to obtain the formula ijSJl, so we still have to extremize 
^ with respect to T and 9 (see discussion in Sec. Ill B|) . 

The classical problem (jA12|) is conveniently dubbed 
T/9 boundary value problem. Eqs. (|A12b|) and HA12c|) 
imply eight real boundary conditions for two complex 
second-order differential equations l|A12a|) . However, one 
of these real conditions is redundant: Eq. l|A12b|) im- 
plies that the (conserved) energy is real, so the condition 
Inipo — > is automatically satisfied (note that the oscil- 
lator energy Eqsc = www = oje^uu* is real). On the other 
hand, the system ljA12p is invariant under time transla- 
tions along the real axis. This invariance is fixed, e.g., 
by demanding that Re(X) takes a prescribed value at a 
prescribed large negative time tg (note that other ways 
may be used instead. In particular, for E < Ei{N) it 
is convenient to impose the constraint ReX{t = 0) =0). 
Together with the latter requirement, we have exactly 
eight real boundary conditions for the system of two com- 
plexified (i.e. four real) second-order equations. 



APPENDIX B: A PROPERTY OF SOLUTIONS 
TO T/e PROBLEM IN FOR THE CASE OF 
OVER-BARRIER TRANSITIONS 

For given E, N there is only one over-barrier classi- 
cal solution which is obtained in the limit e of the 
regularized T/9 procedure. To see what singles out this 
solution, let us analyze the regularized functional 



(Bl) 



where q denotes the variables x{t), x'{t) and T, 9 to- 
gether. The unregularized functional F has a valley of 
extrema q'^{(p) corresponding to different values of the ini- 
tial oscillator phase ip. Clearly, at small e the extremum 
of is close to a point in this valley with the phase 
extremizing Tint [g*" (((C)] . 



(B2) 



Hence, the solution of the regularized T/9 boundary 
value problem tends to the over-barrier classical solution, 
with Tint extremized with respect to the initial oscillator 
phase. 
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Because Ui^itix) > 0, Tint is a positive quantity with 
at least one minimum. In normal situation there is only 
one saddle point of F^, so by solving the T/9 boundary 
value problem one obtains the classical solution with the 
time of interaction minimized. 



APPENDIX C: CLASSICALLY ALLOWED 
TRANSITIONS: ONE-DIMENSIONAL EXAMPLE 

The difhculties with bifurcations of classical solutions 
emerge in quite a general class of quantum mechanical 
models. To illustrate this statement, let us consider the 
case of one-dimensional quantum mechanics, where the 
result is given by well-known WKB formula. We will 
show that the origin of the above difficulties can be seen 
in one-dimensional model also. The implementation of 
the regularization technique is explicit in one dimensional 
case. This makes it easy to see how our technique allows 
one to join smoothly classical solutions relevant to the 
tunneling and allowed transitions. 

In quantum mechanics of one degree of freedom only 
one variable X{t) is present, which describes motion of 
a particle with mass m = I through a potential bar- 
rier U{X). The motion is free in the asymptotic regions 
X — > ±c». The semiclassical calculation of the tunneling 
exponent is performed by solving the classical equation 
of motion 



SS 
SX{t) 



= 



on the contour ABCD in complex time plane, with the 
conditions that the solution is real in the asymptotic past 
(region A), and asymptotic future (region D). The rele- 
vant solutions tend to X —i- — oo and X — > +oo in regions 
A and D, respectively. The auxiliary parameter T is re- 
lated to the energy of incoming state by the requirement 
that the energy of the classical solution equals to E. The 
exponent for the transition probability is 



F = 2Im5 - ET 



(CI) 



One notes the resemblance of these boundary condi- 
tions to the ones on the tunneling coordinate X in two- 
dimensional system. 

In quantum mechanics of one degree of freedom, the 
contour ABCD may be chosen in such a way that the 
points B and C are the turning points of the solution. 
Then the solution is real also at the part BC of the con- 
tour. Indeed, a real solution at the part BC of the con- 
tour oscillates in the upside-down potential, T/2 equals 
to half-period of oscillations, and the points B and C are 
the two different turning points, X = 0. The continua- 
tion of this solution, according to the equation of motion, 
from the point C to the positive real times corresponds to 
the real-time motion, with zero initial velocity, towards 
X +oo; the coordinate X{t) stays real on the part CD 
of the contour. Likewise, the continuation back in time 



from the point B leads to real solution in the part AB of 
the contour. In this way the reality conditions are satis- 
fied at A and D. The only contribution to F comes from 
the Euclidean part of the contour, and one can check that 
the expression IjCljl reduces to 



Xc 

F{E) = 2 J y^2{U{X) - E) dX 



(C2) 



Xe 



which is the standard WKB result. 

The solutions appropriate for the classically forbidden 
and classically allowed transitions apparently belong to 
different branches. As the energy approaches the height 
of the barrier C/q from below, the amplitude of the os- 
cillations in the upside-down potential decreases, while 
the period T tends to a finite value determined by the 
curvature of the potential at its maximum. On the other 
hand, the solutions for E > Uq always run along the real 
time axis, so the parameter T is always zero. Hence, the 
relevant solutions do not merge a,t E = Uq, and T{E) has 
a discontinuity at E — Uq- The regularization technique 
of Sec. IIII Al removes this discontinuity and allows for 
smooth transitions through the point E — Uq. The only 
difference with quantum mechanics of multiple degrees of 
freedom is that in the latter case the bifurcation points 
exist not only at the boundary of the region of classi- 
cally allowed transitions, but also well inside the region 
of classically forbidden transitions (but still a.t E > Es, 
see Introduction and Sec. Ill CI) . 

To illustrate the situation, let us consider an exactly 
solvable model with 

UiX) = 



cosh^ X 



We implement our regularization technique by formally 
changing the potential 



U{X) 



'UiX) , 



(C3) 



which leads to the corresponding change of the classical 
equations of motion. Here e is a real regularization pa- 
rameter, the smallest parameter in the model. At the 
end of the calculations one takes the limit e ^ 0. 

We do not change the boundary conditions in our reg- 
ularized classical problem, i.e., we still require that X{t) 
be real in the asymptotic future on the real time axis 
and that X{t') be real as t' — oo on the part A of the 
contour ABCD. Then the conserved energy is real. The 
sphaleron solution X{t) — has now complex energy (be- 
cause the potential is complex). Hence, the solutions to 
our classical boundary value problem necessarily avoid 
the sphaleron, and one may expect that the solutions 
behave smoothly in energy. 

The general solution to the regularized problem is 



E 



- E 



■ sinh X — — cosh I 
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where to is the integration constant. The value of Im to is 
fixed by the requirement that ImX = at positive time 
t +00, 



iarg[e 



E] 



The residual parameter Re to represents the real-time 
translational invariance present in the problem. The con- 
dition that the coordinate X is real on the initial part AB 
of the contour gives the relation between T and E, 



T 
2" 



1 



{TT + arg E)] 



(C4) 



For e = and E < \, the original unrcgularized result 
r/2 = tt/^/TE is reproduced. 

Let us analyze what happens in the regularized case 
in the vicinity of the would-be special value of energy, 
E ^ Eci = l. It is clear from Eq. lEH) that T IS now a 
smooth function of E. Away from E — \, Eq. ljC4p can 
be written as follows, 



T 
2" 



'2E{E - 1) 



forbidden region, 1 — E ^ e 
allowed region, £^ — 1 e. 

(C5) 



Deep enough in the region of forbidden transitions, when 
1 — -E ^ e, the argument in equation ^VA\ is nearly zero 
and we return to the original tunneling solution. When 
E crosses the region of size of order e around E — 1, the 
argument rapidly changes from 0(e) to — tt, so that T/2 
changes from 7r/V2 to nearly zero. Thus, at > 1 we 
arrive to a solution which is very close to the classical 
over-barrier transition, and the contour is also very close 
to the real axis. This is shown in Fig. O We conclude that 
at small but finite e, the classically allowed and classically 
forbidden transitions merge smoothly. 



At E < 1, the limit e ^ is straightforward. For 
-E > 1 a somewhat more careful analysis of the limit 
e ^ is needed. From Eq. (|C5|I one observes that the 
limit e — > with constant finite T < ■E\pl leads to so- 
lutions with E — 1. Classical over-barrier solutions of 
the original problem with E > Es = 1 are obtained 
in the limit e — > provided that T also tends to zero 
while T = T/e is kept finite. Different energies corre- 
spond to different values of r. And that is what one 
expects — classical over-barrier transitions are described 
by the solutions on the contour with T = 0. 
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